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Abstract 

We construct planar black hole solutions in AdSs and AdS 4 in which the boundary 
CFT is perturbed by marginally relevant quenched disorder. We show that the entropy 
density of the horizon has the scaling temperature dependence s ~ T^d.-i)/z 
d = 2,3). The dynamical critical exponent z is computed numerically and, at weak 
disorder, analytically. These results lend support to the claim that the perturbed CFT 
flows to a disordered quantum critical theory in the IR. 
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1 Introduction 


Consider a quantum field theory in the presence of quenched disorder, that is, spatially 


random couplings. As with simpler spatially uniform couplings, the effects of quenched 
disorder can be relevant or irrelevant mm- When relevant, the disorder can drive the 


quantum theory into various possible nontrivial low energy phases. Of interest to us in 
this paper will be disordered fixed points. At a disordered hxed point, the low energy 
physics exhibits an emergent scale invariance. Such fixed points have a different structure 
to, for example, relativistic conformal field theories because the spatial dependence of the 
couplings at the fixed point means that momentum is not conserved at long wavelengths. 

It is not easy to find controlled instances of disordered fixed points where the disordered 
coupling is stabilized at a finite value [2]. Nonetheless, such critical theories without mo¬ 
mentum conservation are very interesting candidates to understand the universal behavior 
of bad metals [3]. It is therefore of interest to have concrete examples at hand. Evidence 
for a disordered fixed point was recently found in a holographic system [Jj . The disordered 
fixed point itself is dual to a highly inhomogeneous extremal black hole horizon. 

In [1], a GET described by a gravity dual was perturbed by marginally relevant disorder. 
It was shown that the disorder averaged geometry in the far IR regime exhibited an emergent 
scaling invariance 



( 1 ) 


Here ds"^ is the disorder averaged metric. This result was obtained both by resummation 
of perturbation theory and by full-blown numerics. The averaged metric was therefore seen 


to be characterized by a dynamical critical exponent z that determines the relative scaling 


of space and time [5] . However, it was not completely clear what physical quantities would 
be determined by this disorder averaged metric. 


In this paper we will revisit the system considered in [3], but now placed at a nonzero 


temperature. Our main result, that we again obtain both analytically and numerically, is 
that the entropy density of the system scales with temperature as 




( 2 ) 


with d the number of spacetime dimensions of the quantum field theory and with the same 
z appearing as in the disorder averaged zero temperature IR metric 0. Namely, in an 
expansion in small disorder strength V, 


Z = l + l7r‘^/2-ir (I) V^ + O (V^) . 


(3) 
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Thus, we have shown that the averaged metric indeed accurately captures the scaling prop¬ 
erties of a disordered fixed point. In fact, we will see that the temperature scaling of 
the entropy Q is equal to the temperature scaling of the entropy of the averaged metric, 
although the coefficients need not agree. 

The main technical achievements in this paper are the perturbative and numerical con¬ 
struction of the T > 0 disordered black hole spacetimes. ^ describes the general setup. 
^ and ^ obtain the solution perturbatively in disorder and describe a resummation of 
logarithms, ^constructs the solutions numerically. In ^we discuss open questions. 

2 Setup 

In this section we review the holographic description of a CFT perturbed by marginal 
disorder [3]. The starting point is a real scalar coupled to gravity in d-|- 1 dimensions Q The 
action is: 

S = [ d'^+^x ^/=^ [R-A- 2VAd>V^d> - 4y($)] . (4) 

Here = SttGn and A = — is the usual AdS^+i cosmological constant. The resulting 

equations of motion are: 

0 = □$ - H'($), Rab = 2Va^Vb^ + j^9ab [4H(4>) + A] . (5) 

For the scalar potential, we take a negative mass squared: 

4 ^( 4 *)= -^ 4 . 2 . ( 6 ) 

The holographic dictionary (e.g. [H]) tells us that is dual to an operator O with dimension: 



More explicitly, this is seen by considering the asymptotic behavior of the scalar near the 
AdSrf+i boundary. In the Poincare patch, where the line element approaches: 

ds^ = ^ (^rjabdx^'dx’' + dr^ ^ -^ , (8) 

^Throughout our discussion we will let d denote the boundary spacetime dimensions. We will always 
work in signature (—h + • • • -t). For index conventions, capital Latin indices will denote all bulk directions, 
lowercase Latin letters, a, b etc. , will denote boundary spacetime directions, while middle lowercase Latin 
letters, i,j etc. , will denote the boundary spatial directions. 
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as r —)• 0 , the scalar has the following form near the boundary: 

^.(r ^ 0) = + • • • • (9) 


The correspondence then tells us that <hi is identified with the source for O while ^2 encodes 
the response [ 6 ]. 

Our interest in this paper is to consider the effect of a disordered source for O at finite 
temperature T. To be explicit, we will work with a short ranged, quenched, Gaussian 
disorder ensemble, where the ensemble of sources is determined by: 


= 0 , 


- y^). 


( 10 ) 


All other moments of the distribution are then fixed by Wick contraction. Note that, as 
behts quenched disorder, the random sources only depend on the boundary spatial direc¬ 
tions. Our analytic discussion later will involve a resummation of perturbation theory in 
V, whereas the numerics will be exact in V. 

We will focus on the case of ‘marginal’ disorder; that is, we will take the distribution 
to saturate the Harris criterion, which determines when short-range disorder affects critical 
phenomena mm- A simple heuristic way to see this result is to note that since is dual to 
an operator of dimension A, dimensional analysis tells us that [‘hi] = d — A and therefore 


(10) suggests we assign V a dimension of: 

2[$i] = 2[V]+d-l 


[q = ^-A, 


( 11 ) 


We expect then that the disorder is relevant if A < , irrelevant for A > and marginal 

for A = Requiring [H] = 0 fixes the value of /r in ([^ to be: 

d ^-1 


M = 


( 12 ) 


To realize the disorder, we will use a spectral representation [7], writing the source as: 

N-l d-l 

$l(a:*) = V n COs{ki^n,X^ + -fn,) . (13) 

{ni}=l i=l 

Here the 7 ^^ are random phases uniformly distributed over (0, 27r) while the specifics of the 
distribution are determined by the constants and the selection of ki^m- To strictly 

capture the disorder in the thermodynamic limit we must take A —)• 00 . The disorder 
average of a quantity / is then given by: 


/ 


lim 

N^OO 


d-l N-l 

nn 

_i=l ni=l 


2tt 


/. 


(14) 
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We will consider the simplest case of a short range, Gaussian and isotropic disorder distri¬ 
bution, which corresponds to; 


=C=i 2VAk 


d-l 


ki,ni = UiAk , 


Ak-^ 

AT 


(15) 


The wavevectors of the modes making up the disordered source (13) therefore range from 


ko/N to ko- These are the IR and UV cutoffs on the disorder distribution, respectively. In 
principle we could take the spacings Ak to depend on the direction of k, but for simplicity 
we take an isotropic distribution. 

Since we will be working at hnite temperature, it is important to keep the various 
scales in mind. It is useful to consider the two dimensionless parameters: kq = ko/T and 
= kq/N. The spectral representation requires we take iV —)• oo and physically we want 
ko T, but the order of limits is important. Since our aim is to describe a disordered 
system at small but hnite temperature, we should be taking the N ^ oo limit hrst, and so 
in what follows we will work with the following hierarchy: 

ko 


Kir <C 1 <C Ko , i.e. 


N 


< T <C fco . 


(16) 


3 Perturbative geometry 

In this section we perturbatively construct the spacetime deformed by the disordered bound¬ 


ary source (15). This involves solving the bulk scalar held wave equation subject to the 


disordered boundary condition, computing the energy momentum tensor of this scalar held, 
and then computing the backreaction on the metric. Our analytic discussion will largely 
focus on the lowest order thermodynamic corrections. These will be logarithmic in temper¬ 


ature, suggesting a natural resummation. We will show later in section 4.1 that to obtain 
the entropy as a function of temperature to this lowest order, it is sufficient to obtain the 
metric that follows from the backreaction produced by the disorder averaged scalar stress 
tensor. This amounts to hnding the leading disorder averaged correction to the metric, 
which we will now do. 

In the appendix we specialize to the d = 2 dimensional case and determine the correction 
to the geometry for generic scalar conhgurations, without averaging. It will be noted that 
despite the expressions being rather complicated, no essentially new physics is found from 
the full, configuration dependent expressions. 
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3.1 Geometry at 0{V^) 


We will work throughout in the Poincare patch at finite temperature. Therefore in the limit 
P —> 0 , the line element reduces to (from here on we set L = 1 ): 


ds2 


1 


-f{r)dt^ + 


dr^ 

fir) 




2=1 


(17) 


Here /(r) is the emblackening factor /(r) = 1 — (^)'^, where r+ is the horizon radius. In 
terms of the temperature T, we have T = 4 ^^. The entropy density of the thermal state 
in the dual field theory is then given by the familiar Bekenstein-Hawking entropy: 


1 1 


d'^-ix*V=7 


1 1 


1 

AGn 


(Ivr)'^-^ 



(18) 


Here 7 is the metric induced on the horizon from This scaling of the entropy with 

temperature is an important result to keep in mind as our primary objective in this work 
is to determine the modification of this scaling relation due to the disorder. 


3.2 Scalar solutions at OiV) 

We now turn on the disorder with strength V. That is, we introduce bulk scalars whose 


near-boundary behavior, modulo a factor of gives the boundary source (13). This 

scalar solution is determined by the wave equaton in the background ([ItI): 


0 = r‘^+^dr + r^df<A>d) + 


(19) 


Since we are using a spectral representation of the source on the boundary, we decompose 
our bulk scalar into harmonics as well. 


$(i)(^,a.*) = CV'^ (t>kir)Y[cos{ki^n,x'- +'yi) , 
{ni} i 


( 20 ) 


where the (with k = Ife^J) now solve the ODE: 


0 = r 


d+l 


dr [r ^'’fdrfk^ - [k^r^ - n)(t)k 


( 21 ) 


The holographic prescription is to find the linear combination of solutions to this differential 
equation which behave as in ([^ near the boundary and are regular at the horizon . That 
is, due to the various constants that have been factored out, we are to pick the solutions of 
( 21 ) that are regular at the horizon and behave at small r as (/>fc(r —?• 0 ) = + ■ ■ ■ ■ 
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This differential equation does not have a closed form solution for d > 2, as the em- 
blackening factor introduces d singular pointsj^ Fortunately, for reasons that will become 
clear below, we only need the large k behavior of the scalars. These large k modes will 
be responsible for the leading IR singular behavior after disorder averaging. In the large k 
regime we can employ a WKB approximation to find (letting n = kr^ and p = r/r_|_): 




d-1 

P " 


exp 


-i^P^Fi {h\A + h,p' 


( 22 ) 


The WKB limit here is k = kr+ —?• oo, or k/T —?• oo. These modes are largely insensitive 
to the presence of the horizon, decaying well before reaching the horizon, whereas the small 
K <C 1 modes will only weakly vary between the boundary and the horizon. 


3.3 Geometry at 0{V‘^) 


Once the scalars are turned on in the bulk, they source the Einstein equations at order 
As mentioned above, to start with we will find the geometry induced by the averaged stress 
tensor. To leading order this is the disorder-averaged finite temperature metric. To that 
end, we calculate the averaged trace-reversed stress tensor: 




The needed averages are simple to calculate using the spectral decomposition (20), and the 


resulting sources are: 




K%fQ 


rr 




2 '^- 2 (d- 1 ) 



1 / 2(72 

2'=*-2 


E 


f{4>'kr - 


1/2(72 

2"*-2 



P 4>l' 

d — 1 r 2 ’ 


P \ ^ 
d — 1J r2 


(24) 

(25) 

(26) 


Since we have taken an isotropic distribution, the scalar sources in the spatial direction Qa 
are equal for all i. This will simplify the resulting geometry considerably. 

With these sources, we search for a perturbative solution for K <C 1 of the form: 


ds = ^ 


d ^2 


d-1 


f{r) (1 + V^A{r)) dt^ + — + (1 + V^B{r)) 


( 27 ) 


^In the appendix we have the closed form solution for d = 2. 
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Plugging this ansatz into Einstein’s equations then yields the following system of coupled 
ODEs: 


0 = fA" + , 


2r 

0 = + 


2r 

B" - 


r 2 / 


+ 2k,j\tQ- 


N'^rr ) 


0 = /5" - 


, d + f{d-2)^, /./,o..2 


-B - A + 2Kj^Qii. 


(28) 

(29) 

(30) 


r r 

These equations can be re-expressed as a first order constraint equation plus two decoupled 
second order equations. The second order equations are: 




rf.l-d yl/2 


^l-dj3/2 


d+{d-2)f 


dr A 


+ jA{r), 


0 = 2r'^-^fl‘^dr 


n^d - 1 


drB 


+ 3 B{r) 


(31) 

(32) 


where the scalar sources have been repackaged into: 


jA{r) = - 24 I [d + 3(d - 2)f]r^Qtt -[d-{d- 2)f\ (^fQrr - ^ 0n^ I , 

/ + /©rr + ©i 


3B{r) = 


2kj^ 

d-1 


(33) 

(34) 


The decoupled second order equations (31) and (32) can be solved exactly. It is useful 
to work with the rescaled coordinate 


P 


r+ 


(35) 


If we further rescale the sources v'^Ja/b 3a/Bi then we can write the solution as: 
A{p) = ai a2f~^^^[id - 2)f - d] 


+ A 


dp 


d.lp P' 


d-l 


d-{d-2)f{p) d-{d-2)f{p) 


5(p) = /3i + /324 [4'(p)-1 


P/^P) 

+ 1 /■' 


d 


A P' 


id-1 


1 - 


f^np) 

/l/ 2 (p) 


f^^‘^{p)3A{p) 

[d+{d-2)f{p)Y 

3b{p) ■ 


(36) 

(37) 


Thus we have reduced the scalar backreaction to two integrals. There are a number of 


integration constants in (|36|) and (|37|). Indeed, the behavior of the metric at the horizon is 

(or, equivalently, 


entirely determined by these constants. Plugging these solutions into 
into the first order equation that can be derived from the three equations above) shows that 
/32 is fixed by 0 : 2 , /32 oc 0:2 (we only need to say they are proportional, as they will both be 
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zero shortly). In showing that the remaining equation is satisfied, it is important to verify 
that the disorder averaged energy momentum tensor is conserved. 

The constants are fixed by the boundary conditions we impose on the metric, both at 
the conformal boundary and at the horizon. The physical requirements for the geometry 
are that it is asymptotically AdS^+i and that it is regular at the horizon. Regularity at the 
horizon is easily seen to require 0:2 = 0 (and hence (32 = 0). At the conformal boundary, we 
require A(0) = B(0). The actual value of A(0) can be scaled away by redefining coordinates 
so we will impose the simple condition A(0) = B{0) = 0, that is: 


1 

ai = - - 


dp 


r} I nd—l 
^ Jo P L 


d- {d-2)f{p) 


^ 1 = - \ 


dp 


1 - 


f^^Hp) 
1 


- 2 


f^‘^{p)jAip) 

] [d+{d-2)f{p)]^ 


jBip) ■ 


djo p^-n p/^{p)\ 

Thus we obtain an explicit expression for the metric at order V'^. 


(38) 

(39) 


3.4 Large momentum backreaction 

In this section, we will look at the backreaction induced by the large momenta (relative to 


the temperature) scalar modes, where the WKB solutions (22) are valid. In particular, we 
are interested in the behavior of the metric at the horizon, as this is what determines the 
entropy. To get the metric at the horizon, we must determine the integration constants ai 


and /3i (because the integrals in (36) and (37) vanish at the horizon p = 1). We will see 


that these large momentum modes give the leading singular contribution to oi and /3i. 

Since the average stress tensor is a decoupled sum of contributions from each wavevector, 
we write Ja/b = 'T^jk^aib^k- The large momentum contribution to ai,/3i can be readily 
evaluated in the limit k —)■ 00, using the WKB solutions 


0:1, K - 


1 


dp 


d- {d-2)f{p) 


I3i,k - , 


/ 1 

Jo P 


f^Hp) 

^-'^AKr(d +1) 

+ 0{n} 

2M{d- 

1) 

/■' dp 

1 - 

1 

Jo 

p/^{p)\ 


- 2 


pi d^f/'^{p)jA^^{p) 

[d+(d-2)/(p)]2 




JbAp) = o(^k 


(40) 

(41) 


The integrals are performed by noting that at large k, the small p (near boundary) region 


dominates (specifically, p ~ 1/^)- This means that the WKB solution (22) goes like 
while the remaining terms in the integrand can be expanded about p = 0. 

Recall that n = A:r_|_, so that the WKB limit is k/T -A 00. We can sum the above large 
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K terms up to find: 


1^0 p 

fl / 

K>1 


"0 Ki-"r(d + i) 

2^d{d-l) 


f 


dK ■ 


d 

7r2' 


■r{i) 


2k 


~_vrf-ir(01ogKo. (42) 

In the last step we have picked out the singular contribution due to the upper endpoint 


kq S> 1 of the integral. The large k contributions to /3i in (41) are smaller by a power of k 


than the large k contributions to ai. There is no singular contribution in that case. 

Since ai and (5i will also generically receive non-zero contributions from all momenta, 
we write: 


ai = r]i- 7r2 ^T (0 log^o , 


I3i = m- 


(43) 


If r/i, r]2 have no singular dependence on r_|_ in the range kq/N <C 1 <C kq, then this 
information is all we need to determine the low temperature scaling of the entropy density. 
But this last statement is indeed true. It may be verified by numerically evaluating all 
the necessary integrals. Alternatively, we can note physically that the only place that such 
dependence could arise is from the IR cutoff kq/N; however the modes near the IR cutoff, 
in the limit N ^ oo, are essentially constant between the boundary and the horizon and 
so their contribution will be that of a k = 0 mode, which won’t introduce any singular r_|_ 
dependence. 


The logarithmic divergence in (42) has essentially the same origin as the zero tempera¬ 


ture logarithm found in [3], as well as the logarithms arising in earlier works laie]- 

4 Thermodynamics 

We can easily obtain the thermodynamic properties of the averaged geometries constructed 
in the previous section. First we recall that the temperature of the horizon is determined 
by the surface gravity, k, which now receives contributions from ai: 


T = — = ^ 

27r dvrr. 


d 




^ +^f{r+)A{r+)V^ + 0{V^) 

1 + (0 log(/cor+) 


+ 0{V^). 


(44) 

(45) 


If, in the spirit of [1], we throw caution to the wind and exponentiate the logarithm, then 
to this order we can write: 




( 46 ) 
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where 


1 


z = l + -7r2-^r{l)V^+ 0{V^). 


(47) 


This is precisely the dynamical critical exponent identified at T = 0 in [5] . The constant rji 
has gone into the prefactor in this scaling relation. 

Now that we know how the temperature scales with the horizon, we can determine the 
entropy density scaling; 


s = 


1 1 




1 r 


-(d-i) 


AGn V 


d^-^x 


l + ^V^B(r+) + 0(V^) 


(48) 


r=r+ 


Since B{r^) is simply an r_|_-independent constant at low temperatures, we see that this 
entropy scales as: 


-(d-l) 

S ^ T\ ~ i 2 


(49) 


This scaling relation is the first incarnation of our primary result. In particular, we see 
that the exponent z is a true critical exponent and that the disorder has indeed affected 


thermodynamic properties. The result (49) relates the entropy and temperature of the 


averaged metric. In the following subsection we show that this relation also holds for the 
true entropy as a function of temperature. 


4.1 Configuration dependence 

In this section, we discuss the entropy of a typical, configuration dependent metric. Without 
the enhanced symmetry of the averaged configuration to simplify matters, the line element 
will in general look like: 


ds^ = 1 


dr 


/(r) (l + V‘^A{x\r)) dt^ + + ^ (^ + Bij{x\r)) dx'dx^ 


(50) 


We cannot solve analytically for the metric functions A, Bij in general. However, the entropy 
depends only on the induced metric on the horizon. To leading order: 


1 1 


s = 


d^-^x^ = 


4GnV 


1 r 


-{d-i) 


AGn V 


d x\ det(l + V‘^Bij{rj^)) 


1 r. 


r=r+ 

-(d-i) 

+ 


4Gjv V 


d'^-^x 




1 1 
A:Gjv rl“^ 


l + ^V^'^BUr+) 


(51) 
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Here we have used det(l + eA) = 1 + etrA + O(e^) and the fact that the metric components, 
being given by linear combinations of sines and cosines, are self averaging. This result tells 
us that, to the order at which we are working, the entropy density only depends on the 
averaged spatial metric, and we can use the result of the previous section to conclude that 
s ~ just as before. Therefore to determine if the entropy of the averaged metric is 

distinct from the entropy of a typical configuration dependent metric all we need to do is 
find the surface gravity. 

The surface gravity is easily worked out to lowest order in full generality: 

^2 Af ir+)? ^ V^A{r+, X*) + 


K = 


\f{r+)\ 


l + ^V^A{r+,A + 0{V^) 


(52) 


It is a theorem (for metrics of the form we are considering) that the surface gravity must 
be constant along the horizon. Therefore we can replace H(r+,x*) in the previous equation 
by its average: 


K = 




1 + V^A{r+) + 0{V^) 


(53) 


More explicitly, averaging k over the horizon is trivial since it is a constant, whereas aver¬ 
aging H(r_|_,x*) over space is the same as averaging over the disorder ensemble. Using our 
results for the average metric to we can deduce: 


K ~ r ~ 


l/'(^+)l 


1 — ^7r2 ^T (^) (log A;or+-|-const.) 


'f 


+ ’ 


(54) 


where z is again the exponent identified above and in 

The temperature and entropy scalings, r+ ~ and s ~ rA combine to give: 


s r 




(55) 


now as a result for the actual entropy as a function of temperature. In the remainder of the 
paper we will verify this result with full blown numerics, beyond the perturbative regime. 


5 Numerics 

In order to construct the fully backreacted black hole solution, at finite disorder strength 
V, we use the DeTurck trick mm- The method works in a general number of dimensions, 
and we shall use it for d = 2, 3 (boundary spacetime dimensions). We will detail the d = 3 
construction, since it is more involved and has not yet appeared in the literature]^ 

^The DeTurk method was used to construct disordered spacetimes in [4] and m- Other numerical studies 
of strong disorder in holography have been in the probe limit [SIIIIIIS]- 
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The black hole solution we search for is static, which means we can introduce an adapted 
coordinate system in which dt is a Killing direction. Furthermore, the line element should 
be invariant under the discrete transformation t —?■ —t. The most general line element and 
scalar field compatible with these symmetries can be written 


ds^ = ^ 


B 


- {l-y'^)AyX dr + - - .d^T 

I -yi 

y\Si (^dxi + Fidy + F2da;2^ + y\ S2 {dx2 + Fsdy^ 


y+ 


, (56) 


(57) 


where A, B, Si, S 2 , Fi, F 2 , F 3 and <I> comprise a total of 8 functions that depend on y, 
xi and X 2 - The first step in using the DeTurck method is to choose a reference metric g. 
This reference metric should have our desired boundary conditions {i.e. contains a regular 
horizon and has the correct asymptotics). For the reference metric, we choose the planar 


Schwarzschild black hole, which can be obtained from the line element (56) by setting 


A = B = Si = 82 = 1 and Fi = F 2 = F 3 = 0 . The reference metric does not depend on the 
boundary spatial coordinates xi and X 2 , and so is automatically periodic with respect to 
these. Finally, is a parameter which controls the black hole temperature: dvrT = 3?/+. 
The second step in the DeTurck method consists of solving the following set of equations 

1 


□ $ _ y'($) = 0, 


gAB [4K($) + A] = 0 . 


(58) 


where = g^^\^^Q{g) — ^^q{ 9 )\ and r(^) is the Levi-Civita connection associated with 
the reference metric g. Furthermore, K(<h) is given by ([^ with a mass saturating the Harris 
criterion (12). The advantage of this method is that the above Eqs. (58) form a set of 
elliptic PDFs |10j for the metric ansatz (56), unlike the Eqs. ([^. Eor solutions of Eq. (58) 
to also be solutions of Eq. ([^, we must have = 0. In some cases (such as vacuum 
Einstein), there is a proof that all solutions to Eq. (58) also have = 0 [H]. In our case, 
we lack such a proof, so we must verify that = 0 after solving the equations. The local 
uniqueness property of elliptic equations guarantees that solutions with / 0 cannot be 
arbitrarily close to those with = 0 . 

In order to complete the system of partial differential equations, suitable boundary 
conditions must be imposed. In addition, these must be consistent with zero DeTurck 
vector . At the boundary, located at y = 0 , we demand A(0, xi,X 2 ) = i?( 0 ,xi,X 2 ) = 
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5i(0,xi,X2 ) = S'2(0,X1,X2) = 1 and Fi(0,xi,X2) = F2(0,xi,X2) = F5{0,xi,X2) = 0. Fur¬ 
thermore, we demand $( 0 ,xi,X 2 ) = d*i(xi,X 2 ), with the scalar source ‘hi(xi,X 2 ) defined 


in Eq. (13). The reader might be surprised with the extra factor in Eq. (57) dependent 
on 7/_|_. However, we note that, asymptotically, the relation between the Fefferman-Graham 
coordinate r defined in Eq. (|^ and y reads y = yj^r + 0{r'^). 

At the horizon, y = 1, the Einstein-DeTurck equations demand A(l, xi, X 2 ) = H(l, xi, X 2 ), 
which is equivalent to having a well defined bifurcating Killing horizon, with our choice of 
reference metric. The boundary conditions at the horizon for the remaining variables follow 
from expanding the equations in a power series off the horizon - they all turn out to be of 
the Robin type. For the xi and X 2 directions, we demand periodic boundary conditions. 


Now we are in a position to solve the PDE system (58) subject to the above mentioned 
boundary conditions. To solve the equations, we use a standard damped Newton-Raphson 
iteration algorithm based on pseudo-spectral collocation on a Chebyshev (in y) and Fourier 
grids (in xi and X 2 ). In d = 2, there are additional subtleties associated with the boundary 
behaviour of the scalar field <h, which introduces several non-analytic behaviours in the 
metric and scalar field functions. To deal with these, we patch a grid of finite differences 
onto the Chebyshev collocation grid parametrising the holographic radial direction y [Ij. 

The space of solutions is 4-dimensional, depending on V, ko, N and T. However, since 
our underlying UV microscopic theory is conformally invariant, we only need to worry about 
dimensionless ratios of these quantities. In order to access the true IR physics, we need to 


preserve the hierarchy presented in Eq. (16). That is, we must make sure the temperature 
range we probe is in between the short and long distance cutoffs on the disorder distribution. 
In the Schwarzschild background with V = 0, the temperature is given by T = (i/(47rr_|_), 
and it is r_|_ rather than l/T that sets the scale that should be compared to cutoffs. This is 
helpful because it pushes the IR cutoff down to lower temperatures (by a factor of (i/(47r)) 
than the rough window (16) would suggest In all computations detailed in this section 
we measure all quantities in units of /cq (which effectively sets /cq = 1 ) and either we take 
= 50 in d = 2, or iV = 5 in d = 3. T is then allowed to vary freely in the required range 


(16). Note that for these values of N, the Fourier grids must be very dense, having at least 


"^Having the IR cutoff on the disorder be behind the horizon also resolves the following technical issue that 
arises at T = 0. While the disorder is marginally relevant, the homogeneous mode of the scalar is strongly 


relevant. In the energy range (161 the growth of the homogeneous mode is subdominant to the disorder 
physics due to the presence of many higher harmonics. However, below the IR cutoff on the disorder, the 
homogeneous mode will eventually dominate and drive a flow away from AdS. This is an artifact of needing 
to work with an IR cutoff, and can complicate zero temperature numerics, but not the numerics herein. 
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500 points in the periodic direction in d = 2 and 50 in d = Typical profiles for the 
boundary source in d = 2 and d = 3 are depicted in Fig. Q. 


0.2 

0.1 

<ti 

0.0 

- 0.1 

- 0.2 
(a) 

Figure 1: Disordered sources: Plot (a) shows a scalar source as a function of xi ko 
at the boundary. Plot (b) is a density plot of d>i, now in d = 3, as a function of boundary 
directions xi ko and X 2 ko. In both cases we have chosen V = 0.1. The characteristic width 
of the peaks in these plots is determined by the short distance cutoff, Ax|peak ~ T^/ko- 



Having these solutions at hand, there are several quantities we can monitor. We decided 
to focus on the entropy, since it is a direct probe of the infrared geometry. We will discuss 


the results for d = 2 and d = 3 separately, starting with d = 2. In Fig. (2(a)) we show a plot 
of the logarithmic derivative of the entropy, as a function of the black hole temperature, 
for several values of the disorder amplitude V. The data is represented by disks, and the 


solid lines indicate the analytic prediction of Eq. (55), namely s oc with z given 


by (47). From top to bottom, we have V = 0.1, 0.2,..., 1.0. The agreement between our 
perturbative analytic prediction and the numerics is striking. This numerical scaling result 
is compelling evidence for the emergence of a disordered fixed point at T = 0, characterized 
by a dynamical critical exponent z. 

In d = 3 the calculations are more involved, since we have to generate many solutions 
to the 3D PDE system we described above. This means we do not have as much data as 
for the d = 2 case. In particular, we have focussed on a single value V = 0.1. We also have 
a narrower window of temperatures in which to access the IR scaling regime (116^ because 


®The choice of number of points in each periodic direction is such that we should be able to resolve 
up to the fifth harmonic of the highest wave number appearing in our scalar field potential ( |13[ ). Since 
each harmonic descendent decays exponentially m in multiples of the relevant wavenumber, we expect our 
resolution to capture all the relevant physics. 
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N is smaller. In Fig. |2(b)| we plot the logarithmic derivative of the entropy, as a function 
of the black hole temperature for V = 0.1. The disks represent the data, and the solid line 


indicates the analytic prediction s oc of (55), with z again given by (47). Again, 


the agreement as the temperature is lowered is rather encouraging. 
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Figure 2; Emergence of an IR dynamical scaling exponent. Plot (a) shows the 
logarithmic derivative of the entropy for several values of E in d = 2. These plots have 
N = 50. From top to bottom, we have V = 0.1,0.2,..., 1.0. Plot (b) shows the logarithmic 
derivative of the entropy for V = 0.1 in d = 3. This plot has N = 5. 


The computations of the entropy we have discussed are for a given realization of disorder. 
This is the entropy we have been after. With the numerical data at hand, we can compare 
this (physical) entropy with the entropy of the averaged metric, as discussed in previous 
sections. The averaged metric is easily obtained from the numerics by integrating over x 
and y (see the more extended discussion in HD- To compare the averaged entropy with 
the entropy of the averaged metric, we did the following: we computed the entropy of the 
average metric and the entropy of the full metric. We subtracted one from the other and 
found a maximum disagreement around 1%. We then did a similar subtraction, but this 
time for the dynamical critical exponent measured with both entropies, and we found a 
maximum disagreement of 10“^%, which is well within the error in our integration scheme 
in d = 3. We take this as strong evidence that the dynamical critical exponent yields the 
same value whether measured with the entropy of the average or full metric, as we have 
argued in the previous section. This result also substantiates the claim in [1] that the 
averaged metric is a useful bulk quantity to identify the scaling properties of the IR fixed 
point. 
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The reader might also be interested in the difference between the sources shown in 
Figs. Q and the scalar field evaluated at the horizon, For completeness, we show the 
latter in Figs. We see that some of the structure of the UV clearly survives in the IR. At 
first sight, the disorder appears to have been smoothened out on the horizon relative to the 
sources shown in Fig. 0. However, this is simply the fact that upon renormalizing down 
to the horizon, structure on scales smaller than the temperature scale have been integrated 
out. This illustrates the need to keep the long distance cutoff on the disorder distribution 
sufficiently large in order to access the correct disorder physics in the IR geometry. 



Figure 3: Disordered horizons. Plot (a) shows the scalar field evaluated at the 
horizon as a function of xi Uq. Plot (b) is a density plot of <h'^, now in d = 3, as a function 
of boundary directions xi ko and X 2 ko. The sources for these solutions are those shown in 
Fig. Plot (a) is at temperature T/ko = 0.0478 while plot (b) has T/ko = 0.0798. The 
characteristic width of the peaks in these plots is now determined by the temperature scale, 
so that Ax I peak ~ vrr+. This is the expected statement that the temperature serves as the 
short distance cutoff on the disorder distribution at the horizon. 


Finally, in an attempt to characterize the geometry of the horizon, we plot in Fig. @ 
the Ricci scalar, of the induced metric on a spatial cross section of the horizon - a 

disordered horizon. 


6 Discussion 

In this paper we have presented evidence for the existence of a disordered fixed point in the 
far IR of a spacetime with a marginally relevant disordered boundary source. In addition 
to constructing numerical disordered black hole spacetimes, we showed that the dynami- 
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Figure 4: Disordered horizon - Ricci scalar. Plot of the Ricci scalar of the induced 
metric on a spatial cross section of the horizon. The parameters used are the same as in 
Fig-i Because the metric depends on the square of the scalar field, the metric functions 
oscillate twice as quickly and hence the structures appear half the size of those in Fig. 

cal critical exponent z of the IR fixed point could be obtained by resumming logarithmic 
divergences that appear in perturbation theory in the strength of the disorder. However, 
this perturbation theory is an expansion about the UV spacetime. The whole point of fixed 
points is that they are self-contained and well-defined without reference to a UV comple¬ 
tion. Therefore, an intrinsic description of the disordered horizon (at T = 0) as a solution 
to Einstein’s equations should exist. Characterizing the disordered horizon on its own, IR, 
terms could potentially lead to a greatly expanded understanding of what extremal hori¬ 
zons can look like. It would, presumably, explain why a naive resummation of logarithmic 
divergences at low orders in perturbation theory appears to give the correct answer for the 
dynamical critical exponent. It would also clarify exactly which quantities can be accurately 
determined from the corresponding disorder averaged spacetime. 

Given the construction of the background geometries, it is very natural to study correla¬ 
tion functions in these backgrounds. If the full frequency and momentum dependence can be 
found, then this should verify the value of z that we have found, giving correlators that are 
scaling functions of ujjk^. The study of transport, in particular, in these backgrounds was 
one of the motivations to construct these solutions in the first place. The disordered fixed 
point does not conserve momentum and will therefore have intrinsically finite transport 
properties [3]. This is a qualitatively distinct regime from the case in which disorder can 
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be described as an irrelevant perturbation about a clean IR fixed point [nillBlIIHlEolEI]. 

Finally, having developed the numerical and analytical methods needed to understand 
disordered spacetimes, we may soon be in a position to tackle the more difficult case of 
relevant (rather than marginally relevant) disorder. 
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A d = 2 solution 

In d = 2, one can go (slightly) further in describing the geometry of individual configurations 
in terms of some explicit indefinite integrals. The hrst key difference is that in d = 2 the 
scalar solution at finite temperature can be written explicitly: 

N-l 

$(l)(r,x) = (j)n{r) cos{knX + Jn) , (59) 

n=l 

where: 

Ur) = C^Ilr (I - (;^) . (60) 

Here Pu{x) is the Legendre function (for the particular values of the index v = —^ + ik 
these are also known as parabolic cylinder functions). Throughout this section we will use 
the notation: 


= knX + 7n , 


^nm — • 


(61) 


A further simplification in the d = 2 case is simply that the spatial metric is much less 
complicated when there is only one spatial direction. The order geometry is given as: 


ds^ = ^ 


dr^ 

-/(r) (l -F V‘^A{r, x)) dd -h + (l + E^R(r, x)) dx^ 

J[r) ^ 


(62) 
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The equations of motion for A, B are particular combinations of Einstein’s equations: 


0 = rf-^/‘^dr 

-1 

CO 

to 

2 r 
H-2 


r 



0 = 2/^/2 

0 = r^f-y^dr —drA +r‘^dlA-rfdrA-rdrB-A^x{^^^^f. 


0 = rf^'^dr 

cq 

cd 

1_ 

2 r 

H-2 


r 



j ,/2 


Plugging in the scalar solution (59) and (60), it is straightforward to determine: 


A{r, x) = ai{r) + / ^/‘^{r)a 2 {x) + A{r, x ), 


(63) 

(64) 

(65) 

( 66 ) 

(67) 


where: 


A{r,x) = al 


m 


V)cos 29m + 


E' 

m^n 


iir)kn 


cos 




6 t,„ ^ cos 6 ^ 


hr 


( 68 ) 


where the radial profiles for the harmonics are: 


alnir) = f f dh , (69) 

J r 

Omn = 2/"^/^(r) [ df f/^{f)^m{h)(j)n{r) ■ (70) 

Jr 


The function ai{r) in (67) is given by integral expressions similar to what we found for 


general d: 


0-1 (r) =r]i+ r] 2 r 1 ^/^(r) - 1 

rr+ 


T3 1 -\/f{r)/f{r) h‘^f{h){(i>mir)) - (b^A:^ +/i)(^^ . (71) 


Turning now to the spatial metric component, it can be written as some integration 
functions and a particular solution: 


B{r,x) =/3i{x) + fd 2 {x)r‘l. /^/^(r) - 1 +B{r,x). 


(72) 


The particular solutions can be written as: 


B(r, x)=Y^ [b'^^ir) + 6)„(r) cos^ 9^] + Y1 [^mn(^) cos 0+„ + bl,^{r) cos 9^ 

m m^n 


(73) 
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where: 


C(r) = 2 

+ 

1 

/l/2(^) _ 

J^l/2 j'j:' 


[fkm(l)m{r)f , 


(74) 

blnir) = 2 

+ 

1 

[/l/2(^) _ 
J^l/2 

- 


'klr + tl)(lL{f)] , 

(75) 

bLnir) = 

+ 

1 

[/l/2(^) _ 
J^l/2 



kmkn4^m{x) 4^n{x)~\ 






-h(t)m{r)4>n{r)] , 


(76) 

blinir) = 

+ 

1 

[/l/2(^) _ 
^*1/2 


{r^ [f{f)(t>'^{h)(i)'n{r) 

+ kmkn4^m{x) (pn{x)~\ 



-^J-(t>ra{f)(|)n{r)] • (77) 


There are a few integration functions left we need to fix. In particular, the constants 
r/i, ry2 and the functions a 2 {x), I3i{x), f32{x). As we discussed before, r]i receives a logarithmic 
contribution from the WKB modes, which works out to be; 

m = - ^y‘^^ogkor+ H-, (78) 

while r] 2 , 02 and /32 are all forced to vanish by regularity at the horizon. This leaves /3i(x), 
which is determined by fixing A(0,x) = B{0,x). Our solutions for A and B allow one to 
write explicit integral expressions for /3i, though they do not appear to be illuminating so 
we have not written them out. 
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